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Abstract. The filamentation instability (FI) of counter-propagating beams of 
electrons is modelled with a particle-in-cell simulation in one spatial dimension 
and with a high statistical plasma representation. The simulation direction is 
orthogonal to the beam velocity vector. Both electron beams have initially equal 
densities, temperatures and moduli of their nonrelativistic mean velocities. The FI 
is electromagnetic in this case. A previous study of a small filament demonstrated, 
that the magnetic pressure gradient force (MPGF) results in a nonlinearly driven 
electrostatic field. The probably small contribution of the thermal pressure gradient to 
the force balance implied, that the electrostatic field performed undamped oscillations 
around a background electric field. Here we consider larger filaments, which reach a 
stronger electrostatic potential when they saturate. The electron heating is enhanced 
and electrostatic electron phase space holes form. The competition of several smaller 
filaments, which grow simultaneously with the large filament, also perturbs the balance 
between the electrostatic and magnetic fields. The oscillations are damped but the final 
electric field amplitude is still determined by the MPGF. 
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1. Introduction 

The filamentation instability (FI) driven by counterpropagating electron beams amplifies 
magnetic fields in astrophysical and solar flare plasmas [1-5] and it is also relevant for 
inertial confinement fusion (ICF) [6] and laser-plasma interactions in general [3 Ej. It 
has been modelled with particle-in-cell (PIC) and Vlasov codes [9-17] taking sometimes 
into account the ion response and a guiding magnetic field. It turns out that the FI 
is important, when the beam speeds are at least mildly relativistic and if the beams 
have a similar density [18] . Otherwise its linear growth rate decreases below those of 
the competing two-stream instability or mixed mode instability [19] . 

The saturation of the FI is attributed to magnetic trapping [9] . More recently, it has 
been pointed out [T3J, [H] that the electric fields are also important in this context. An 
electric field component along the beam velocity vector Vf, is driven by the FI through 
the displacement current. This component is typically weak and its relevance to the 
plasma dynamics is negligible compared to that of the magnetic and the electrostatic 
fields. The FI is partially electrostatic during its linear growth phase, if the electron 
beams are asymmetric due to different densities. Symmetric electron beams result in 
purely electromagnetic waves with wavevectors k _L vt, [T9J, [201. A nonlinear growth 
mechanism is provided in this case by the current of the electrons, which have been 
accelerated by the magnetic pressure gradient force (MPGF). 

The electromagnetic and electrostatic components separate in a ID simulation box, 
because the gradients along two directions vanish in the Maxwell's equations. The 
electrostatic field is polarized in the simulation direction, while the electromagnetic 
components are polarized orthogonal to it. If both electron beams have an equal density 
and temperature, the electrostatic field component along the wavevector k can only 
be driven nonlinearly. We select here a direction of our ID PIC simulation box that 
is orthogonal to vj,, through which this nonlinear mechanism can be examined in an 
isolated form. The equally dense and warm counterstreaming beams of electrons have 
the velocity modulus \vj,\ = 0.3c. The ions are immobile and compensate the electron 
charge. The mildly relativistic relative streaming speed ~ 0.55c implies, that the growth 
rate of the FI is significant. At the same time, any relativistic mass changes can be 
neglected during the growth phase and the saturation of the FI. 

The initial conditions of the plasma equal those in the Refs. [21| [22]. The size 
distribution of the filaments could be sampled with the help of the long ID simulation 
box in Ref. [21]. A pair of current filaments, which are small according to this size 
distribution, has been isolated in Ref. [22] . It could be shown that the electrostatic field 
is indeed driven by the MPGF for this filament pair. The electrostatic field performed 
undamped oscillations around a background one. The latter excerted the same force on 
the electrons as the MPGF. Here we assess the influence of the filament size. 

This paper is structured as follows. Section 2 discusses briefly the PIC code, the 
initial conditions and the key nonlinear processes. The results are presented in the 
section 3, which can be summarized as follows. The electrons are heated up along the 
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wavevector k by their interaction with the wave fields. As we increase the filament size 
the peak amplitudes grow, which are reached by the magnetic and by the electrostatic 
field when the FI saturates. The electron heating increases with the filament size and 
large electron phase space holes form, which interact with the electromagnetic fields of 
the filamentation modes. The large box sizes allow the growth of more than one wave 
and the filamentation modes compete. The electrostatic field oscillations are damped 
or inhibited and the amplitude modulus converges to one, which equals that expected 
from the MPGF. We confirm that the strength of the electrostatic force on an electron 
is comparable to that of the magnetic force, when the FI saturates. The extraordinary 
modes are pumped by the FI |14| . The results are discussed in section 4. 

2. The PIC simulation, the initial conditions and the nonlinear terms 

The PIC simulation method is detailed in Ref . [23] . Our code is based on the numerical 
scheme proposed by [21]. The phase space fluid is approximated by an ensemble of 
computational particles (CPs) with a mass m cp and charge q cp that can differ from 
those of the represented physical particles. The charge-to-mass ratio must be preserved 
though. The Maxwell-Lorentz equations are solved. The plasma frequency of each beam 
with the density n e that we model is tu p = (e 2 n e /m e e )°' 5 and Q, p = y/2u p . The electric 
and magnetic fields are normalized to E^ = eE/cm e Q p and B^ = eB/m e Q p . The 
current is normalized to Jn = J/2n e ec and the charge to pn = p/2n e e. The physical 
position, the time and speed are normalized as x^ = x/X s with X s = c/fl p , t N = tfl p 
and vn = v/c. The normalized frequency ojn = oj/VL p . We drop the indices N and 
x, t, u, E, B, J and p are specified in normalized units. The equations are 



with p cp = m cp T cp v cp . Here v cPiX is the component along x of v cp . The currents 
jcp oc q C pV C p of each CP are interpolated to the grid. The summation over all CPs gives 
«J, which is defined on the grid. The J updates E and B through (1). Our numerical 
scheme fulfills (2) as constraints. The new fields are interpolated to the position of each 
CP and advance its position x cp and p cp through (3). All components of p are resolved. 

Two spatially uniform beams of electrons with q cp /m cp = —e/m e move along z. 
Beam 1 has the mean speed u&i = i>& and the beam 2 has = — i>&i with u 6 = 0.3. 
Both beams have a Maxwellian velocity distribution in their respective rest frame with 
a thermal speed v t h = ^{k^T '/m e )°' of Vb/v t h = 18. The negative electron charge is 
compensated by an immobile positive charge background. The initial conditions are 
p, J, E, B = 0. Figure [H displays the k spectrum of the unstable waves. The growth 
rates of the FI modes are close to the maximum value, while relativistic effects are 
still negligible. The growth rate spectrum with ku = relevant for our simulations 
peaks with 5m = 0.29 at kj^K ~ 10. A filamentation mode with /cmA s = 7 has been 



V x E = -d t B , V x B = J + d t E, 

V -E = p, V • B = 0, 

dfPcp = q cp (E[x cp \ + v cp x B[x cp \) , d x x cp = v, 



(1) 
(2) 
(3) 
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Figure 1. (Colour online) The growth rates in units of fl p as a function of the 
wavenumber in the full k space, where A s fcii (A s fcj_) points along (orthogonal) to Vb- 
The growth rates of the FI modes with h\ = are comparable to that of the two-stream 
mode with kj_ = and to those of the oblique modes. The growth rates for kit = 
decrease to zero for k± — > and they are stabilized at high k± by thermal effects. The 
growth rate maximum for ku = is 5m — 0.29 and it is reached at fejwA s ~ 10. 

considered in detail previously [22], while we investigate here larger filaments. The box 
length Li = 2 for the simulation 1 and the filamentation mode with ki = 2n/Li grows 
at the exponential rate 0.92 5m- The box length of the simulation 2 is L 2 = 2.8 and the 
growth rate of the filamentation mode with k 2 = 2tt/L 2 is 0.86 5m- The growth rates 
decrease rapidly for lower k and these modes are no longer observed in PIC simulations 
[2 lj . Both simulations resolve x by N g = 500 grid cells with the length A x and use 
periodic boundary conditions. The phase space distributions f\(x,v) of beam 1 and 
f 2 (x,v) of beam 2 are each sampled by N p = 6.05 ■ 10 7 CPs. The total phase space 
density is defined as f(x, v) = fi(x, v) + f 2 (x, v). 

Each electron beam constitutes prior to the saturation of the FI a fluid with 
the index j, which has the density n,j(x) = f v fj(x,v)dv and the mean velocity 
Vj(x) = f v vfj(x, v)dv. The normalized momentum equation for such a fluid is 

d t {n jVj ) + VirijVjVj) = -VPj - n 3 E + V(BB) - VB 2 / 2 + B x 9 t E, (4) 

where the thermal pressure tensor Pj is normalized to 2m e n e c 2 . The restriction to one 
spatial dimension implies, that the gradients along y and z vanish. The FI results in 
this case in the initial growth of B y and of a weaker electric E z . The thermal pressure 
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is initially diagonal due to the spatially uniform single-Maxwellian velocity distribution. 
The x-component of the simplified fluid momentum equation is 

dt(njVj,x) + d x (rijV? x ) = -v 2 th d x n j - n j E x - B y d x B y + B y d t E z . (5) 

The thermal pressure gradient v 2 h d x rij is valid, as long as the electron beams have 
not been heated up. Let us assume that the displacement current and the thermal 
pressure gradient can be neglected, leaving us with the term njE x and the MPGF as 
the key nonlinear terms. The fluid momentum equations can be summed over both 
beams and we consider the right hand side of (J5j). As long as E x is small, the electron 
density is not spatially modulated and n\ + n 2 pa 1. The nonlinear terms cancel out, if 
E x = —2Byd x B y . It could be demonstrated for a short filament in Ref. [22] that this 
is the case, even when the FI just saturated. The E x oscillated in time and after the 
saturation with the amplitude E B = —B y d x B y around a time-stationary E B . 

3. Simulation results 

3.1. The scaling of B y , E x and E B with the box length 

The beam velocity Vj, \\ z and the electrons of both beams and their micro- currents are 
re-distributed by the FI only along x. The initially charge- and current-neutral plasma 
is transformed into one with J z (x,t) ^ 0. The gradients along the y, ^-direction vanish 
in our ID geometry. Ampere's law simplifies to d x B y = J z + d t E z , resulting in the 
growth of B y and E z . The MPGF drives E x . The B x = in the ID geometry and 
E y , B z remain at noise levels. The right-hand side of ([5]) depends on E x , E z and B y , as 
well as on their spatial gradients, which should vary with the filament size. 

We want to gain qualitative insight into the scaling of the field amplitudes with the 
filament size and determine if E x is driven by the MPGF also for the large filaments. 
The fields that grow in simulation 1 and 2 are compared to those discussed previously 
in Ref. [22J that used the box size L c = 0.89. Figure [2] shows the respective dominant 
Fourier component of B y , of E x and of 2E B . The amplitude moduli of the mode with 
k s = 2tt/L s are considered for B y and those of the 2k s mode for E x and 2E B . The 
subscript s is 1, 2 or c and refers to the respective simulation. The amplitudes of 
By increase with an increasing box size. After the FI has saturated, we find that 
B y (ki,t) ~ 2B y (k c ,t) and B y (k 2 ,t) pa 2.5B y (k c , t). The increase of the saturation value 
of B y (k s ,t) with L s is consistent with magnetic trapping [9]. The magnetic bouncing 
frequency lui, = (vi,k s B[k s , t}) 1 ^ 2 in our normalization. The FI should saturate once Ub is 
comparable to the linear growth rate of the FI, which is approximately constant for the 
box sizes L c , L\ and L 2 (Fig. [I]). A lower k s supports a larger B y (k s ,t). The Uf, ~ 0.2 
for simulation 1 is comparable to the linear growth rate Ui pa 0.25. 

After the saturation, the E x (2k 1)2 ,t) > 2E x (2k c ,t) and E x {2k h t) > E x (2k 2 ,t). The 
E x (2ki,t) > 3E x (2k c ,t) while Li/L c pa 2.2. The electrostatic potential in simulation 1 
is thus larger by a factor 6, which should result in a more violent electron acceleration 
than in the box with the length L c . The thermal pressure gradient force is potentially 
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Figure 2. (Colour online) Panel (a) compares the B y (k s ,t) and panel (b) the 
E x (2k s ,t) in the boxes with the size L%, Li and L c (dashed curve). £ , s(2fc s , t) (dashed 
curve) is compared with E x (2k s ,t) (solid curves) for the box size L\ (c) and Li (d). 



more important for larger filaments and it may modify the balance between the 
nonlinearly driven E x and the MPGF. However, an excellent match between E x (2ki,t) 
and 2Es{2ki,t) is observed for t < 50, due to which the two nonlinear terms on the 
right hand side of (jSJ) practically cancel for simulation 1. The E x (2k2,t) ~ 2EB(2k2,t) 
in simulation 2 for 30 < t < 42 and for 46 < t < 53. Both fields disagree in between 
these time intervals and a local minimum is observed. The field and electron dynamics 
is now examined in more detail for the box lengths L\ and L2. 

3.2. Simulation 1: Box length L\ = 2 

Figure [3] displays the evolution of the relevant field components. The B y (x,t) rapidly 
grows and saturates at t ~ 45. It is initially stationary in space but it oscillates in 
time until t pa 65, which implies that B y (x,t) does not immediately go into its stable 
saturated state. The B y (x,t) shows only one spatial oscillation and the filamentation 
mode with the wavelength k\ = 2tc/Li is thus strongest. However, the interval with 
the large positive B y (x,t pa 45) covers < x < 0.9, while that with the large negative 
B y (x, t pa 45) is limited to 1.2 < x < 1.7. This mode is thus initially not monochromatic. 
The saturated structure formed by B y (x,t) drifts after t ~ 65 to lower x at a speed 
< 0.01 and it remains stationary in its moving rest frame. The E z (x,t) grows initially 
in unison with B y (x,t) and it is shifted in space by 90° with respect to B y (x,t), as 
expected from Ampere's law. Oscillations of E z (x,t) are spatially correlated with those 
of the B y (x,t) for 45 < t < 65. The E z (x,t) undergoes a mode conversion at t pa 65 
into a time-oscillatory and spatially uniform E z (x,t). Figure [3](c) demonstrates that 
E x (x,t) is following the drift of B y (x,t) towards decreasing x, but that its wavenumber 
is twice that of B y (x,t). The B y (x,t) is stationary in its moving rest frame, while 
E x (x,t > 70) is oscillating around an equilibrium electric field with an amplitude and 
spatial distribution that resembles E B (x,t) in Fig. [3]^d). The electric and the magnetic 
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Figure 3. (Colour online) The field amplitudes in the box L\. The panels (a-d) show 
By, E Zl E x and Eb, respectively. The amplitude of B y reaches a time-stationary 
distribution, which convects to decreasing i at a speed < 0.01. The E z and E x 
components are oscillatory in space and in time. The E z is phase-shifted by 90° 
relative to B y when the fields saturate at t » 45. The E x and the Eb are co- moving 
and E x oscillates in time around a mean amplitude comparable to Eb for t > 70. 
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Figure 4. The relevant part of the amplitude spectrum E x (k, t) is displayed for low k 
in (a) and (b) shows that of Es(k,t). The wavenumbers are expressed in units of k\. 
The amplitude moduli of the dominant modes are displayed for k — 1k\ in (c) and its 
first harmonic with k = 4fci in (d), where the dashed curves correspond to Eb- 



forces are comparable in their strength, but their distribution differs. 

Figure H] compares in more detail the moduli of the amplitude spectra E x (k,t) and 
E B {k,t). The amplitudes of the strongest modes fulfill E x (2k\,t) ~ 2E B {2ki,t) until 
t = 50 (See also Fig. [2]). The E x (2kx,t) thus overshoots £ , s(2A;i,t) and it oscillates 
around it after t = 50. The oscillation is damped and the amplitudes of E x (2k\,t) and 
EB(2ki,i) converge. The full spectra E x (k, t) and Esi^k, t) reveal that the mode k = Ak\ 
is also important for t > 100. It probably is a harmonic of the mode with k = 2k± and 
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Figure 5. (Colour online) The 10-logarithm of the phase space densities in units of 
CPs at the time t = 50 (a-c) and t = 120 (d-f) in the box L±: Panels (a,d) show the 
total phase space density f(x,p z ) with the beam momentum po — m e VbT(vb). The 
phase space density fi{x,p x ) of beam 1 is shown in (b,e) and the f2(x,p x ) of beam 
2 in (c.f). Both beams are spatially separated and (e,f) reveal cool electron clouds 
immersed in a hot electron background with momenta of up to ~ po- 



not an independently growing FI mode. Otherwise we would expect that the mode with 
k w 3k\ also grows. The amplitude of E x (4k\,t) is close to that of Es{4:ki,t). 

A dissipation mechanism for the interplaying J x and E x is present, which causes 
the damping and the convergence of E x (x,t) to E B (x,t). The damping persists after 
t = 65, when B y is quasi-stationary in its moving reference frame. The term B y dtE z in 
(EJ) could, in principle, be one dissipation mechanism. However, even at t m 50 when 
dfE z is largest and B y has developed in full, this term is weaker by more than one order 
of magnitude than the MPGF and the term rijE x in simulation 1 (not shown). If the 
term B y dE z would be the damping mechanism, this should have resulted in a noticable 
field damping also in the short simulation box with length L c . A damping of E x (x,t) 
has not been observed in Ref. [22]. The thermal pressure gradient force may provide 
this damping and we examine now the electron phase space density distribution. 

Figure [5] displays the phase space distributions fi(x,p x ) of the beam 1 and the 
f2{x,p x ) of the beam 2 at the times t = 50 and t = 120. The total phase space density 
f(x,p z ) is shown for the same times. The beams reveal a high degree of symmetry 
already at t = 50 and the filament centres are shifted along x by Li/2. The phase space 
structures in Fig. E(b,c) are, however, different at the filament boundaries x ~ 0.5 and 
x ~ 1.5. This difference is responsible for the deviation of the initial B y (x,t) from a 
sine curve in Fig. [3]^a). The phase space distribution at late times reveals, that the 
electrons are heated along p x but not along p z . The filament drift to lower x is visible 
from Figs. E(a,d) and agrees with the observed one of B y (x,t) in Fig. [3j The electrons 
are accelerated along x to a peak speed ~ i>&, which is more than twice that observed 
in the box with the length L c [22J. The peak electron kinetic energy due to the velocity 
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Figure 6. The field amplitudes in the box L^: The panels (a-d) show B y , E z , E x 
and Eb = —B y d x B y , respectively. The amplitude of B y reaches a steady state value, 
which convects to increasing x at a speed < 0.01. The E z and E x components are 
oscillatory in space and in time. The E z is phase-shifted by 90° relative to B y when 
the fields saturate at t « 50. The E x and the Eb are co- moving and E x (x) oscillates 
in time around a mean amplitude comparable to Eb for t > 100. 



component along x thus increases by a factor, which is comparable to the increase in 
the electrostatic potential as we go from a box with length L c to one with L\. This 
strong electron heating is likely to result in higher thermal pressure gradient forces. 
The expression & x n\{x) J v x f(x,v x )dv x has been evaluated (not shown) at t = 75 and 
its peaks reach values ~ 0.1, which are comparable to the MPGF. The width of these 
peaks is small compared to the electron skin depth. 

Movie 1 animates in time the 10-logarithmic phase space distributions fi{x,p x ) 
and fi(x,p z ) of the beam 1 in the simulation 1. The formation of the filaments is 
demonstrated. We observe a dense untrapped electron component immersed in an 
electron cloud that has been heated along the simulation direction by the saturation 
of the FI. The spatial width of the plasmon containing the dense bulk of the confined 
electrons in fi(x,p x ) oscillates in time. The overlap of the filaments in Fig. [5](e,f) is 
thus time dependent and related through its current J x (x,t) to the oscillating E x (x,t) 
in Fig. [3](b). The phase space distribution fi(x,p x ) reveals small-scale structures (phase 
space holes) that gyrate around the centre of the filament. These coherent structures 
result in jumps in the thermal pressure. 

3.3. Simulation 2: Box length L 2 = 2.8 

Figure [6] displays the fields that grow in the simulation with the box length L 2 = 2.8. The 
growth rate map in Fig. [1] demonstrates that the FI can drive simultaneously several 
modes in the simulation box. The mode with k 2 = 2tt/L 2 ~ 2.25 has, for example, 
a lower growth rate than that with k w 2k 2 . We observe consequently oscillations 
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Figure 7. The relevant part of the amplitude spectrum E x (k, t) is displayed for low k 
in (a) and (b) shows that of The wavenumbers are expressed in units of k^. 

The amplitude moduli of the dominant modes are displayed for k — 2&2 in (c) and its 
first harmonic with k = 4fc2 in (d) , where the dashed curves correspond to Eb ■ 



in B y (x,t) along x, which are a superposition of several waves with a k > k 2 during 
the initial growth phase 40 < t < 50. These oscillations merge and only one spatial 
oscillation of B y (x,t) and, thus, a single pair of filaments survive after the saturation 
at t ~ 50. The magnetic field structure convects to increasing values of x at a speed 
< 0.01, but it is stationary in its rest frame after t « 70. The phase of E z (x,t) is 
shifted by 90° with respect to B y (x,t) for 40 < t < 60. The oscillations of E z (x,t) 
undergo a mode conversion during 60 < t < 100 and we observe undamped oscillations 
with k — for t > 100. The amplitude of these oscillations is modulated on a long 
timescale. The E x (x,t) and the B y (x,t) show no correlation until t ~ 70. Thereafter 
the spatial amplitude of E x (x,t) oscillates in time around E B (x,t). The force on an 
electron imposed by E x (x,t) is comparable to that imposed by VbB y (x,t). 

A more accurate comparison of E x (x,t) and E B (x,t) is again provided by the 
moduli of their spatial amplitude (Fourier) spectra, E x (k,t) and E B {k,t). Figure [7J 
displays E x (k, t) and E B (k, t) and compares in more detail E x (2k 2 , t) with E B (2k 2 , t) as 
well as E x (4k 2 ,t) with E B (4k 2 ,t). The amplitudes E x (2k 2 ,t) « 2E B (2k 2 ,t) during the 
exponential growth phase of the FI for 25 < t < 45 (See Fig. [2]), the amplitude moduli 
then have a local minimum and continue to grow after this time. We identify the likely 
reason from E B (k,t) in Fig. [7(b). The E B (3k 2 ,t) competes with E B (2k 2 ,t) at t ~ 50. 

A large amplitude modulus of E B (3k 2 ,t) evidences that B y (x,t) is not a sine wave 
at this time. If B y oc sin(fc 2 x), then E B oc sin {k 2 x) cos (k 2 x) and E B (k,t) would be 
composed of a wave with k = 2k 2 . The periodic boundary conditions would also allow 
for a By oc sin (2k 2 , t) and here E B would involve a wave with k = 4k 2 . An E B (3k 2 , t) can 
thus not be connected to a single filamentation mode. During the linear growth phase 
of the FI prior to t ~ 40, the J z (x, t) can form structures with a wideband wavenumber 
spectrum (See Fig. [[]) and their associated B y can grow independently. 
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Figure 8. (Colour online) A time-interval of E z (x, t) and the 10-logarithm of its power 
spectrum PEz{k,uj) are displayed in (a) and (b). Wavenumbers are given in units of 
fe. Peak 1 is at lu < 0.5 and k — &2- Peak 2 is observed at k = k-i and w»l and peak 
3 at k — and uj w 1. The Byikz^t) is shown in (c), the E z (k2,t) in (d) and E z (0,t) 
in (e), all normalized to the maximum of B y (k2,t). 



Once the MPGF in Eq. [5] has reached a significant strength, the FI saturates. 
The strength of the MPGF increases with k, due to the larger d x B y (x,t) of the rapid 
oscillations. The B y oc sin (k 2 x) should maximize the magnetic field strength for a given 
MPGF. This may explain why this mode is the dominant one after t = 70 despite its 
lower growth rate. The decrease of EB(2k 2 ,t) in Fig. W(d) at t ~ 45 is tied to the 
saturation of EB{3k 2 ,t). The E x (k,t) in Fig. [7(a) has a broadband spectrum within 
50 < t < 75, which is probably caused by the current J x arising from the rearrangement 
of the filaments. After this time, E x (2k 2} t) « E B (2k 2 ,t) and E x (Ak 2 ,t) « E B (4:k 2 ,t). 
The E x (2k 2 ,t) does not show oscillations around Es{2k 2l t) as the simulation 1. The 
filament rearrangement inhibits an oscillatory equilibrium between J x and 

Figure [8] examines the mode conversion of the electromagnetic E z component 
observed in Fig. 0(b). The PEz{k,uj) is the squared modulus of the Fourier transform 
of E z (x,t) over space and over 45 < t < 100. The dispersion relation shows three 
peaks. Peak 1 has a k = k 2 and u < 0.5 and it is tied to the E z (x,t) of the FI mode. 
This mode grows exponentially and aperiodically. Its frequency spectrum is thus spread 
out along <jj. Its energy can leak into the peak 2 at k = k 2 and lu ~ 1. The E z (x,t) 
is orthogonal to B y (x,t) and peak 2 corresponds to an extraordinary mode, similar 
to the slow extraordinary mode. Peak 3 has a k = and uj ~ 1 and it corresponds 
to a spatially uniform oscillation in an extraordinary mode branch. The intermittent 
behaviour of E z (x,t) in Fig. El(a) results in a broadband spectrum in k and uj. These 
turbulent wave fields can couple energy directly to the high-frequency electromagnetic 
modes and excite a discrete spectrum if the boundary conditions are periodic [14] . 

The interplay of the waves belonging to the three peaks in Fig. [8](a) is assessed 
with the moduli of the amplitude spectra B y (k 2 ,t), E z (k 2 ,t) and E z (0, t) in Figs. [S)(c- 
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Position Position 

Figure 9. The 10-logarithmic phase space densities in units of CPs at t = 50 (a-c) 
and t = 120 (d-f) in the box L2: Panels (a,d) show the total distribution f(x,p z ) 
with po = rn e VbT(vb). The beam temperature along p z is unchanged. The distribution 
fi(x,p x ) of beam 1 is shown in (b,e) and the f2(%,Px) of beam 2 in (c,f). The electrons 
of both beams spatially separate and (e,f ) reveal a dense electron component immersed 
in a tenuous hot electron background, which reaches a thermal width ss pq. 



e). The B y (k2,t) and E z (k2,t) grow at the same exponential rate until they saturate 
at t « 50, evidencing that they belong to the same FI mode. The B y (k2,t) maintains 
its amplitude after t = 50, while E z (k2,t) decreases until t ~ 120 and remains constant 
thereafter. The E z (0,t) grows in the same time interval to its peak amplitude, which 
suggests a parametric interaction between these modes. The amplitude modulation 
in Fig. [8(e) must be caused by a beat between two waves, which are similar to the 
slow- and fast extraordinary modes in the limit k = 0. Both modes are undamped on 
the resolved timescales. One may interpret the parametric interaction as a three-wave 
coupling between the waves corresponding to the peaks 1-3 in Fig. [8](b), resembling 
the system of Ref. [25J. However, here the B y (x,t) varies spatially and the parametric 
interaction may involve more of the waves of the spectrum in Fig. E(b). 

Figure [9] displays the phase space densities fi t 2(x,p x ) and f(x,p z ) at the times 
t = 50 and t = 120. Figure [9](a) demonstrates that the electrons of both beams have 
been rearranged by the FI. The filaments have not yet reached the stable symmetric 
configuration, because the most pronounced density minima at x ~ 1.5 for beam 2 and at 
x ~ 2.5 for beam 1 are not shifted by L 2 /2. This asymmetry results in the Eb{^2, t) 7^ 
and in the broadband E x (k,t) at this time in Fig. The spatial gradients of B y (x,t) 
and E x (x, t) are high at t ~ 50 and the Lorentz force changes rapidly with x, explaining 
the complex phase space structuring in Fig. E3(b,c). The phase shift of L 2 /2 of the 
density maxima of both beams has been reached at t = 120 in Fig. [9(d). The B y (x,t) 
is stationary in its rest frame at this time in Fig. E(a). The electrons are heated up 
from an initial thermal spread of p x /po ~ 0.05 with po = ^eVb^{vb) to a peak value of 
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Px ~ Po i n Fig- Efof). The mean momentum of each beam varies along p z in response 
to a drift imposed by E x (x,t) and B y (x,t) but no heating is observed in this direction. 

Movie 2 shows the 10-logarithmic phase space density projections fi(x,p x ,t) and 
fi(x,p z ,t) of beam 1 in the simulation 2. It demonstrates that only the core electrons in 
Fig. [9]remain spatially confined. The heated electrons, which have in some cases reached 
a momentum p x that is comparable to the initial beam momentum, are untrapped. The 
heated electrons move practically freely and they ensure that the beam confinement is 
not perfect. The trapped electrons maintain the J z (x, t) ^ and, thus, the B y (x, t) ^ 0. 
The trapped electrons slowly move to larger values of x. The associated shift of J z (x, t) 
causes the slow drift of B y (x,t) in Fig. E(a). The movie visualizes the formation of the 
phase space beams and their evolution into phase space holes in f\{x,p x ). 

4. Discussion 

We have examined here the electron beam filamentation instability (FI) in one dimension 
and in an initially unmagnetized plasma with immobile ions. The FI has been driven 
by nonrelativistic symmetric electron beams with the same initial conditions as those 
considered previously [211 E2] . The electric field along the one-dimensional box, which 
is oriented orthogonally to the beam velocity vector, can only be generated nonlinearly 
if the beams are symmetric [20]. The fluid equations show that the relevant nonlinear 
mechanisms can be the magnetic pressure gradient force (MPGF), the thermal pressure 
gradient force and a term due to the displacement current. The magnetic tension may 
become important in multi-dimensional simulations, but not for initial conditions similar 
to ours [26]. The term due to the displacement current is weak in our simulations. 

It has been observed in Ref. [22] that the electrostatic field performs after 
the saturation of the FI undamped oscillations around a time-stationary background 
electric field. The amplitude of the oscillatory and of the background electric field 
are both given by EB(x,t) ~ —B y d x B y . The phases of both fields are fixed such, 
that E x (x,t ) = 2E B (x,t ) at the saturation time t . This amplitude ensures that 
the nonlinear terms due to the MPGF and due to the electrostatic field cancel each 
other approximately in the fluid equations when the FI saturates. The thermal pressure 
gradient force did not visibly contribute in the simulation of the small filament pair |26j, 
possibly because of the only modest heating of the initially cool beams. Here we have 
assessed the importance of the filament size with the help of two ID PIC simulations, 
which used two different box lengths that were larger than that of the ID box in Ref. 
[221 [26]. The initial conditions for the plasma were otherwise identical. 

We summarize our findings as follows. We have demonstrated for both simulations, 
that E x (x,t < to) ~ 2Es(x,t < t ) during the full exponential growth phase and not 
just at the saturation time t . The FI thus adjusts the electrostatic field during its 
exponential growth phase such, that the dominant nonlinear terms cancel each other. 
Magnetic trapping states that the FI saturates, when the magnetic bouncing frequency 
is comparable to the linear growth rate. The exponential growth rates for the two 
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simulations considered here and that in Ref. [22J are close. The amplitude reached by 
the magnetic field prior to its satuation thus increases with the box length. We found 
that the electrostatic potential driven by the MPGF is 5-6 times stronger for the box 
sizes used here than for the short box in Ref. |22j, while the initial mean kinetic energy 
of the electrons is the same. Consequently, the electron heating is stronger and the 
plasma processes more violent for large filaments. Magnetic trapping is, however, not 
the exclusive saturation mechanism. The electrostatic forces are comparable in strength 
to the magnetic forces when the FI saturates [T3l [14"] . 

The electrostatic field during the intermittent phase has differed in our two 
simulations from that observed in Ref. [22] • The movies demonstrated that this phase 
involves the formation of large nonlinear structures (phase space holes) in the electron 
distribution, which can result in steep gradients of the thermal pressure and in the 
generation of solitary (bipolar) electrostatic wave structures that are independent of the 
fields produced by the FI. The thermal pressure gradient force is comparable to that 
of the other nonlinear terms, but only over limited spatial intervals. The electric field 
component along the beam velocity vector has undergone a mode conversion. Its energy 
leaked into the high-frequency electromagnetic modes [T4"] . 

The wavenumber spectrum of the electrostatic field correlated well with that of 
the MPGF in simulation 1, but the peak electric field overshot the expected one. 
The electrostatic field performed damped oscillations around Eb and both converged 
eventually to the same value. The wavenumber spectrum of the electrostatic field in 
simulation 2 deviated from that of the MPGF in the intermittent phase. Its wavenumber 
spectrum was broadband, while that of the MPGF was quasi-monochromatic. The 
amplitude modulus of the electrostatic field at the wavenumber, which corresponds to 
the dominant Fourier component of the MPGF, jumped to the value expected from the 
MPGF. It did not overshoot and it was non-oscillatory. 

Both simulations here have evidenced that the magnetic field driven by the FI 
organized itself such, that we obtained one oscillation in the simulation box after the 
intermittent phase. This is remarkable, because the exponential growth rate of the 
fundamental wavenumber is below that of its first harmonic. Long waves excert a lower 
MPGF for a given amplitude and the dominance of the fundamental wavenumber may 
thus result from the lower nonlinear damping of this mode compared to that of its 
harmonics. The mode with the fundamental wavenumber considered in Ref. [22] has 
a higher growth rate than its harmonics and the absent mode competition may have 
facilitated the undamped oscillations around the equilibrium. However, the amplitude 
of the electrostatic field in the two simulations discussed here eventually converged to 
that expected from the MPGF and Eb is thus a robust estimate for the electrostatic 
field driven by the MPGF for the considered case. This robustness explains, why a 
connection between the electrostatic field and the MPGF has been observed in a 2D 
PIC simulation [26], where no equilibrium can be reached due to the filament mergers. 

This estimate does, however, not apply if positrons are present. Their current 
reduces that of the electrons. If equal amounts of electrons and positrons are present, 
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the electrostatic field driven by the MPGF is suppressed alltogether [27]. Mobile protons 
will react in particular to the stationary electric field [14J and they will modify through 
their charge modulation the balance between the electrostatic field and the MPGF. 
Highly relativistic beam velocities will probably also modify the balance between the 
MPGF and the electron currents it drives. We leave relativistic beams to future work. 
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